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Abstract 

We study the pattern of activated trajectories in a double well system without 
detailed balance, in the weak noise limit. The pattern may contain cusps and 
other singular features, which are similar to the caustics of geometrical optics. 
Their presence is reflected in the quasipotential of the system, much as phase 
transitions are reflected in the free energy of a thermodynamic system. By 
tuning system parameters, a cusp may be made to coincide with the saddle 
point. Such an anomalous cusp is analogous to a nonclassical critical point. 
We derive a scaling law, and nonpolynomial 'equations of state', that govern 
its bifurcation into conventional cusps. 

PACS numbers: 05.4fl.-a, 05.70.Fh, 42.15.Dp 
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The optimal trajectory concept has been widely used in the theory of noise-activated 
transitions 0-0. In the weak noise limit, when transitions between stable states become 
exponentially rare, one or at most a few trajectories in the system state space are singled out 
as escape paths of least resistance. Also, between the energetically lowest stable state and 
any other state there are at most a few dominant activated trajectories. Such optimal 
trajectories, which are determined by a 'least energy expended' or 'least action' variational 
principle, are experimentally observable In systems that have the property of detailed 
balance, they are time-reversed relaxational trajectories. But in nonequilibrium systems, 
which lack detailed balance, the optimal trajectory pattern extending from a stable state 
may be more complicated. 

Optimal trajectories are similar in many ways to the rays of geometrical optics, which 
characterize in the short wavelength limit the waves emanating from a point source. That is 
because optical rays may be computed variationally too, from a 'least optical depth' princi- 
ple. In a medium with inhomogeneous index of refraction, it is common for a ray family to 
bounce off a caustic surface, leaving the region behind in shadow: not illuminated, or illu- 
minated only indirectly. Other singular features with a catastrophe-theoretic interpretation 
may be produced |||9|]. In a noise-driven system in which detailed balance is violated, the 
pattern of optimal trajectories may contain similar features []I||j,rU|. See Fig. [I]. 



To understand the crossing of optimal trajectories, it is useful to look at the quasipotential 
of the noise-driven system. If e is the noise strength (e.g., e = kT in thermal systems), and 
p(x) denotes the stationary probability density of the system at state x, then a quasipotential 
W = W(x) may be defined phenomenologically by 

p(x) ~ const x exp[— W(x)/e], e — > 0. (1) 

This definition makes sense whether or not the system dynamics are conservative, and 
whether or not the noise acts so as to preserve detailed balance. W equals zero at the 
energetically lowest stable state(s), and IF(x) is essentially the minimum energy needed 
to excite the system to state x. It may be computed as a line integral along the optimal 
trajectory extending to x. 

Formally, W is multivalued at any state, such as the states near a caustic, that is reached 
by more than one optimal trajectory emanating from the energetically lowest state(s). But 
by (0), the least value is dominant, and the trajectories giving rise to others are unphysical. 
The state space of a noise-driven system is typically partitioned by 'switching surfaces', on 
which dominance switches between branches of W . 



FIG. 1. How a caustic (solid curve) could be formed as the envelope of a family of outgoing 
optimal trajectories. 
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The switching of dominance resembles a first-order phase transition in a condensed mat- 
ter system. The similarity is unsurprising, since phase transitions with classical critical 
exponents also have a catastrophe-theoretic interpretation M. Consider, for example, a fer- 
romagnetic system with extensive order parameter m (magnetization), in a magnetic field h. 
Its thermodynamics are determined by a free energy function \I/ = ^(T, h). Below the critical 
temperature T c , \l/ and m = d^fr/dh are multivalued. If the phase transition is classical, i.e., 
of mean-field form, ^/ is three- valued in a sharp-tipped region of the (T, h) plane bounded 
by 'spinodals' of the form h ~ ±const x (T c — T) 3 / 2 . That is because the leading terms in 
the Legendre transform ^W(T, m) = hm — \I/(T, h) are of Ginzburg-Landau type: 

¥ h \T,m) wCj(T- T c )m 2 /2 + C 4 m 4 /4. (2) 

In the catastrophe-theoretic sense ||||, the spinodals are fold caustics. Each is the projection 
of a fold in the graph of which is a two-dimensional surface, onto the (T, h) plane. The 
critical point (T c , 0) from which the spinodals extend is a cusp catastrophe: the projection 
of the point on the graph of \I/ at which the two folds join. A first-order phase transition 
line, on which dominance switches between branches of extends from (0,0) to (T c ,0). 

Switching lines in two-dimensional noise-driven systems are clearly analogous to first- 
order phase transition lines, and caustics to spinodals. Caustics typically terminate at cusps, 
and switching lines also frequently terminate at cusps. So cusps, which are very common, 
are analogous to second-order critical points 0. They are physically important because at 
any cusp, the prefactor 'const' in ([!]), which in general is x-dependent, diverges. 

In previous work we pointed out that in many noise-driven two-dimensional double 
well systems without detailed balance, a cusp may be moved to coincide with the saddle 
point between the wells, by tuning parameters. If they coincide, the Kramers (e — > 0) limit 
of noise-induced interwell transitions is greatly affected. The prefactor in the Kramers 
transition rate formula becomes anomalous: it acquires a negative power of e. 

Precisely at criticality, we were able to approximate the quasipotential W near any such 
'anomalous cusp'. Our expression differed from the polynomial 'normal forms' of conven- 
tional catastrophe theory. In thermodynamics, it would correspond to a nonclassical phase 
transition. In catastrophe theory, it would be interpreted as a nongeneric catastrophe: one 
of the few such of physical relevance to have been discovered since the work of Berry and 
Mount on the short wavelength limit of scattering ||. 

In this Letter, we extend Ref. f|| by analysing the 'unfolding' of an anomalous cusp in 
a typical two-dimensional noise-activated system, as a parameter is moved toward or away 
from criticality. We explain how it may bifurcate into conventional cusps. Our scaling law for 
the bifurcation yields a corresponding law for the divergence of the Kramers prefactor JTT] . 



Consider the following double well model, which is similar to models of blocking dynam- 
ics in glassy systems, where a particle is coupled to a randomly fluctuating barrier whose 
position is coupled to the particle motion |L^]. Let x (a particle position variable) and y 
(a barrier state variable) be nontrivially coupled in such a way that the values ±1 for x and 
for y are stable. If x and y are overdamped and are driven by white noise of strength e, 
their joint dynamics could be modeled by Langevin equations 



x — A„ 



x 



(l - x 2 ) - axy 2 ] + e 1/2 7] x (t) 



y = -\\ y \(l + x 2 )y + e l /\(t). (3) 
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FIG. 2. The pattern of optimal trajectories near the saddle, in the 'broken phase' (a > a c ). 
Here X x = 1, |AJ = 1.1, a = 4.64, and the critical value a c equals 4.62. In the shaded re- 
gion the quasipotential W is multivalued, and a switching line extends from (0, 0) to the cusp 
(x c ,0) m (0.048,0). 

The parameters > and X y < determine the time scales on which x and y evolve, and 
govern the all-important relaxational behavior near the saddle point (0,0), where (x,y) ~ 
(X x x, —\\y\y). The forcing terms (r) x ,r) y ) are a pair of independent Gaussian white noises, 
so that (T)i(s)r]j(t)) equals 5ij5(s — t). The parameter a controls the absence of detailed 
balance: only when a equals yu = | A y | / is there detailed balance, since only in that case is 
the drift field derived from a potential. 

Our results are insensitive to the details of the coupling between x and y, so long as the 
model is symmetric through x = and y = 0. To compute the pattern of optimal trajectories 
emanating from the bottom of the x < well or the x > well, we use the fact that in any 
multidimensional noise-driven system with vector Langevin equation x = u(x) + e 1 / 2 ^^), 
the optimal trajectories are really zero-energy Hamiltonian trajectories, generated by the 
Wentzell-Freidlin Hamiltonian [13| 

tf(x,p) = p 2 /2 + u(x)-p. (4) 

That is because the associated Hamilton's principle is 

5 J L(x, x) dt = 5 J |x — u(x)| 2 dt = 0, 

which is clearly a 'least energy expended' principle. The conjugate momentum p = dLjd~x. 
equals x — u, which measures the system's motion against the drift. 

Figure was obtained from (|3|) and (|j) by integrating Hamilton's equations outward, 
at zero energy, from (1, 0), i.e., from the bottom of the right-hand well. A small portion of the 
x < half-plane is reached by optimal trajectories, but the rest is in shadow. In phase space, 
which is four-dimensional, the optimal trajectories trace out a two-dimensional manifold, 
called a Lagrangian manifold. This manifold lies 'above' only a small part of the x < 
half-plane. It folds over, covering the shaded portion of the (x, y) plane more than once. 
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In the shaded region, the momentum p and the quasipotential W, which equals / p • dx, are 
two- valued (x < 0) or three- valued (x > 0). At any point x, p equals VW(x). 

In Fig. |2|, the parameter a is chosen to be slightly greater than a certain critical value, a c . 
If there is detailed balance, the optimal trajectory pattern contains no singular features, but 
if a is increased through a c , a cusp (x c , 0) emerges from the saddle point at (0, 0) and 
moves toward (1,0). This phenomenon is not peculiar to the model defined by In any 
symmetric two-dimensional double well system that violates detailed balance and has a 
tunable parameter, a similar cusp may be born. The focusing of optimal trajectories at the 
cusp resembles the focusing of rays in a radially symmetric optical system. 

The value a c can be computed from the second-order variational equation 5 2 J L dt = 0, 
which is a criterion for bifurcation. On physical grounds, when a < a c , 5 2 J Ldt computed 
along the on-axis trajectory to the saddle is positive, but when a > a c , it is negative. In 
the model (§), « c turns out (cf. Ref. @) to equal 2/i(/i + 1). 

The cusp (x c , 0) that is present when a > a c resembles a second-order critical point. W is 
three-valued in the sharp-tipped region extending from it, which is bounded by 'spinodals' 
of the form y w ±const x [x c — x) 3 ^ 2 . (See Fig. ^|.) Moreover, there is a switching line 
extending from the saddle at (0, 0) to (x c , 0). As noted, this line is analogous to a first-order 
phase transition line. 

What remains to be understood is how the cusp is born at a = a c . In a three-dimensional 
space with coordinates (x,y;a), there is a line of second-order critical points in the y = 
plane that extends from (0,0; a c ) to (l,0;+oo). By analogy with thermodynamics, one 
might expect (0, 0; a c ) to be a third-order critical point. At any fixed a > a c , the leading 
terms in the Legendre transform W^ y \x,p y ) = yp y — W, close to the cusp, are known to be 
of Ginzburg-Landau type ||[5) : 

W (y) (x,p y ) « C 2 (a)[x - x c (a)}p 2 y /2 + C 4 (a)p 4 y /4. (5) 

One might expect that the correct three-dimensional generalization would be a higher-degree 
polynomial in x, p y , and a — a c . That would allow the birth of the cusp to be viewed as a 
classical phase transition, or one of the generic (polynomial) elementary catastrophes 

In Ref. f2j we presented initial evidence against this. At criticality (a = a c ), we were 
able to construct a scaling solution for W, valid near the x-axis close to the saddle. The 
equation satisfied by the scaling function contained non- integer powers: in fact, powers that 
depended continuously on the model parameter /i. 

By linearizing Hamiltonian dynamics near the saddle, we have now characterized fully 
the behavior of the quasipotential W near a singular point like (0,0;a c ). Our chief new 
result is a cubic equation satisfied by the double Legendre transform of W. It defines a 
higher-order, but nonclassical, critical point. We have also extended our a = a c scaling law 
to the case when \a — a c \ is nonzero but small. These results should extend to any symmetric 
double well system with a tunable parameter. 

Figure^, which was obtained in the 'unbroken phase' (a < a c ), sheds light on behavior 
near criticality. The crucial feature is the two caustics, relics of which appeared in Fig. |2|. 
They form part of the boundary of the 'illuminated' region. Each caustic extends from a 
cusp, which is located very close to the y-axis separatrix between the two wells. As a — > a c ~, 
the cusps neck down to the saddle. Any further increase in a causes the x-axis cusp to be 
born, and to move toward positive x. 
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FIG. 3. The pattern of optimal trajectories near the saddle, in the 'unbroken phase' (a < a c ). 
Here A^ = 1, \X y \ = 1.1, a = 4.61, and the critical value a c equals 4.62. As a — > a c ~ , the two 
cusps close to the y-axis neck down to the saddle. 

The merged cusp at (0,0), when a = a c , is truly anomalous. It is the projection of a 
point on the boundary of the Lagrangian manifold, rather than of a point in its interior. So 
it is a boundary catastrophe: a singular point of a sensitive kind. To explain its bifurcation 
into cusps on the x-axis or y-axis, we must approximate the optimal trajectory pattern in a 
neighborhood of the saddle, for \a — a c \ nonzero but small. 

We accordingly linearize Hamilton's equations, which any optimal trajectory must satisfy, 
near the point (x, y; p x , p y ) = (0, 0; 0, 0) in phase space. A simple analysis (cf. Ref. ||) shows 
that this fixed point has two stable directions, e s = (0, 1; 0, 0) and e s = (1, 0; — 2X X , 0), and 
two unstable directions, e u = (1,0; 0,0) and e u = (0, 1; 0, 2|A y |). The zero-momentum 
directions (no tilde) are eigendirections for relaxational trajectories, which follow the drift. 
In the linear approximation, any optimal trajectory near (x, y) = (0, 0) must satisfy 



where the C's are trajectory- specific constants. 

We now index the 'fan' of optimal trajectories that approach the saddle point, as in Figs. 
H] and [3], by s. The normalization of this index variable is somewhat arbitrary. A reasonable 
choice would be for it to denote distance from the x-axis (at a fixed x > 0, near the saddle). 
With this choice, s = will correspond to the uphill optimal trajectory that climbs toward 
(0, 0) along the positive x-axis. If each coefficient in (||) can be expanded in s about s — 0, 
then by symmetry considerations 



(x,y; Px ,p y ) ss C s e |A ^e s + C e Xx % 
+ C u e Xxt e u + C u e^% 



(6) 



C s = aqs + a 3 s 3 + • ■ ■ , 
C s = b + b 2 s 2 + ■■■, 
C u = C2S 2 + C4S 4 + • • • , 
C u = dts + d 3 s 3 H . 



(7b) 
(7c) 
(7d) 
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We identify the passage through criticality, as a is increased through a c , with the passing 
through zero of the coefficients c<i and d\. So, setting 5 = a c — a, we take c<i and di to be 
linearly proportional to S, to leading order. 

Eq. (Q) comprises four scalar equations. Eliminating s and t among them, we can derive 
'equations of state' relating 5 and any three of the phase space coordinates x, y, p x , and p y . 
When \5\ <C 1, the equation relating x, y, p y , and 5 (near the x-axis), and the equation 
relating x, y, p x , and 5 (near the ?/-axis), turn out to be, respectively, 

= ( Py -2\X y \y) 3 (8a) 

+ hdx^ipy - 2\Xy\y) + kax^py 
= (p x + 2X x x) 3 (8b) 

+ £i8p 2 f- 2 y 2 (p x + 2X x x) + ioyVf 3 



Here fci, fco, and £ are positive constants. At criticality (5 = 0), (|Sa|)-(|8"b"|) reduce to the 
equations we previously obtained by an altogether different technique JJ. 

By definition, the cusp (x c , 0) is the point on the positive x-axis where W or (p x ,p y ) = 
VW stops being multivalued as a function of (x,y), as x increases from 0. It is easy to 
verify from fl5a|) that when 5 is small and negative (i.e., a — a c is small and positive), the 
cusp location x c satisfies x c oc (— 5) l l 2 ^. The 5-dependence of the parent cusps, which have 
y = ±y c , can be computed from (j3b| ). They are the points close to the y-axis where (p x ,p y ) 
first becomes multivalued, as \y\ increases from zero. When 5 — > + (i.e., a — > a c ~), we find 
that the parent cusps neck down at the rate y c oc 5 3 ^ 2 ~ fl . That is so when 1 < \i < 3/2, 
at least; for other /a, the prediction is that there are no parent cusps, and no necking down. 



All these predictions have been numerically confirmed |TT . 

Despite the continuously varying exponents, the emergence of the x-axis cusp is sur- 
prisingly similar to a second-order phase transition. Recall that close to the ferromagnetic 
critical point defined by (^), scaled magnetization M and scaled magnetic field H are related 
by 

M 3 ± M - H = 0, (9) 
or equivalently by the scaling law M = 4>±(H). Here M = Am/\T — T^ 1 / 2 and H = 



Bh/\T - T c | 3 / 2 , with A = ^jC A /C 2 and B = A 3 /C A . The plus (minus) applies when T - T c 
is positive (negative). Eq. ( jSa]) may be rewritten in the form (||), provided that one defines 

M = {p y - 2\\ y \y)/\k 1 5x 2 » + k x^\ 1/2 (10) 
H = 2\X y \k x^y/\k 1 5x 2 ^ + k x^f /2 . (11) 

The plus (minus) applies when k\8x 2 ^ + k^x^ is positive (negative). The law M = <f>±{H), 
in which 5 appears implicitly, provides a unified description of the x-axis behavior both at 
criticality (5 = 0) and away (5^0). 

The most striking consequence of this approach is a general scaling law, showing how 
the quasipotential varies as (x,y;a) moves away from (0,0;a c ), in any direction. It can 
be written using the double Legendre transform W^ x,y \p x ,p y ) = x • p — W, which equals 
/ x • dp. Near the saddle, 



7 



W{x, y) « W(0, 0) - X x x 2 + \X y \y 2 (12) 
W^\p x ,p y ) * -W(0,0) -p 2 x /4\ x +pl/4\\ y \. (13) 

Let i? = R(p x ,p y ) denote the difference between W^ x ' v ^ and the right-hand side of (PT5|) . 
Since W^*) = / x • p dt, R may be expressed in terms of s and t by employing @ and (^). 
By eliminating s and t from the formulae for R(s, t), p x (s, t), and p y (s, t), we find to leading 
order in 5 

R 3 + m x 5p^p 2 y R - m Q pfp y = 0, (14) 

where mi and tuq are positive constants. This is the extension to 5 ^ of a formula derived 
in Ref. §. 

In phase transition language, the cubic equation (|H]) fully characterizes the nonclassical 



structure of the critical point (x, y; a) = (0, 0; a c ), i.e., (p x ,P y ', «) = (0, 0; a c ). Equation (|I4D 
is also interesting from a catastrophe theory point of view. Nongeneric catastrophes, which 
are difficult to classify, may in general be perturbed in an infinite number of ways so as to 
yield generic catastrophes M. But ( |H| ) describes the bifurcation of a nongeneric catastrophe 
(the anomalous cusp) into conventional cusps, in a unique, physically determined way. 

This research was supported in part by NSF grant PHY-9800979. A portion was com- 
pleted while the authors were in residence at the Aspen Center for Physics. 
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